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Abstract 

We introduce from an experimental point of view the main concepts of fluctuation theorems for work, 
heat and entropy production in out of equilibrium systems. We will discuss the important difference between 
the applications of these concepts to stochastic systems and to a second class of systems (chaotic systems) 
where the fluctuations are induced either by chaotic flows or by fluctuating driving forces. We will mainly 
analyze the stochastic systems using the measurements performed in two experiments : a) a harmonic 
oscillator driven out of equilibrium by an external force b) a colloidal particle trapped in a time dependent 
double well potential. We will rapidly describe some consequences of fluctuation theorems and some useful 
applications to the analysis of experimental data. As an example the case of a molecular motor will be 
analyzed in some details. Finally we will discuss the problems related to the applications of fluctuation 
theorems to chaotic systems. 
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1 Introduction 

This article is a review of the main experimental applications of Fluctuation Theorems (FTs) and summarizes 
the plenary talk given at STATPHYS24. In order to define the main contents let us consider several simple 
examples. The simplest and most basic out of equilibrium system is a thermal conductor whose extremities are 
connected to two heat baths at different temperatures, as sketched in fig. 1. The second law of thermodynamics 
imposes that in average the heat flows from the hot to the cold reservoir (from H to C in fig. 1). However the 
second law does not say anything about fluctuations and in principle one can observe for a short time a heat 
current in the opposite direction. What is the probability of observing these rare events ? As a general rule 
when the size of the system decreases the role of fluctuations increases. Thus from an experimental point of 
view it is reasonable to think that such rare events can be observed in systems that are small. A good candidate 
could be for example the thermal conduction in a nanotube whose extremities are connected to two heat baths 
[1], exactly in the spirit of fig. 1. In reality in this kind of experiments the measure of the mean quantities [1] 
is already difficult and of course the analysis of fluctuations is even more complicated. However there is an 
electrical analogy, shown in fig. lb), of the thermal model of fig. la). Let us consider an electrical conductor 
connected to a potential difference V = Va — Vb and kept at temperature T by a heat bath . If the mean 
current I = V/R(R being the electrical resistance of the conductor) is of the order of 1CP 13 A and the injected 
power is about ~ 10~ 19 J (ks is the Boltzmann constant) then the instantaneous current inside the 

resistance has fluctuations whose amplitude is comparable to the mean, as shown in fig. lc). The variance of 
these fluctuations is 5I 2 ~ JcbT/(R tq) where To is the characteristic time constant of the electrical circuit. 
In the specific case of Fig. lc) the current reverses with respect to the mean value. The probability of having 
those negative currents have been studied both theoretically and experimentally in ref. [2, 3] within the context 
of fluctuation theorems, that we will present in sec.3. 

We discuss a second example where the source of fluctuation is not the coupling with a thermal bath, as in 
the case of the electrical conductor, but it is either a chaotic flows or a chaotic force produced by the non-linear 
interaction of many degrees of freedom of a dissipative system sustained by an external driving. We will refer 
to them as chaotic systems. Let us consider a turbulent wind flowing around an object as sketched in fig.2a),b). 
The wind exerts a mean force Fq on the object but the instantaneous force, plotted in fig.2c), is a strongly 
fluctuating quantity which presents negative values [4], i.e. the object moves against the wind,fig.2b). In such a 
case the mean work done by the wind on the object is about 0. J~ 10 20 &;£T and obviously thermal fluctuations 
do not play any role but so does the chaotic flow, which produces the fluctuations. Other similar examples can 
be found for example in shaken granular media [5, 6], discussed in sec. 6. 

These examples stress that in the two experiments the electrical conductor and the turbulent wind we may 
observe the counterintuitive effect that the instantaneous response of the system is opposite to the mean value, 
in other words the system has an instantaneous negative entropy production rate. This effect is induced by 
the thermal fluctuations in the first case and by the chaotic flow in the second case. The question that we 
want analyze in this article is whether the Fluctuations Theorem (defined in section 3) is able to predict the 
probability of these rare events in both cases, i.e. for the stochastic and the chaotic systems. We will take an 
experimentalist approach and we will use experimental results in order to introduce the main concepts. 
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Figure 1 : a) Schematic representation of a conductor whose extremities are in contact with two heat baths at 
temperature Tjj and Tq with Th > Tq- b) Electrical analogy. A conductor of electrical resistance R and kept 
at a temperature T is submitted to a potential difference V = V a — V&. c) Instantaneous current / flowing into 
the resistance using R = 10 MO, T = 300 K and r = 2 ms. 




displacement 



Figure 2: a)and b) Schematic representation of an object suspended by an elastic beam and submitted to the 
pressure of a turbulent wind a) average behavior b) rare event, c) Time evolution of the measured instantaneous 
force exerted by the turbulent wind on the object. The details of this experiment can be found [4] 

The largest part of the article concerns stochastic systems described by a Langevin dynamics. For chaotic 
systems we will mainly discuss the difficulty of comparing the experimental results with the theoretical pre- 
dictions. The article is organized as follows. In section 2 we present the experimental results on the energy 
fluctuations measured in a harmonic oscillator driven out of equilibrium by an external force. In section 3 
the experimental results on the harmonic oscillator are used to introduce the property of Fluctuation Theorems 
(FTs). In section 4 the non linear case of a Brownian particle confined in a time dependent double well potential 
is presented. In section 5 we introduce the applications of the FT, and as a more specific example we describe 
the measure of the torque of a molecular motor. Finally in section 6 we discuss the chaotic systems and we 
conclude in section 7. 

2 Work and heat fluctuations in the harmonic oscillator 

The choice of discussing the dynamics of the harmonic oscillator is dictated by the fact that it is relevant for 
many practical applications such as the measure of the elasticity of nanotubes[7], the dynamics of the tip of an 
AFM [8], the MEMS and the thermal rheometer that we developed several years ago to study the rheology of 
complex fluids [9, 10]. 
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Figure 3: a) The torsion pendulum, b) The magnetostatic forcing, c) Picture of the pendulum, d) Cell where 
the pendulum is installed. 

2.1 The experimental set-up 

This device is a very sensitive torsion pendulum as sketched in fig. 3a). It is composed by a brass wire (length 
10 mm, width 0.5 mm, thickness 50 fim) and a glass mirror with a golden surface, fig. 3c). The mirror (length 
2 mm, width 8 mm, thickness 1 mm) is glued in the middle of the brass wire. The elastic torsional stiffness 
of the wire is C = 4.65 ■ 10" 4 N.m.rad -1 . It is enclosed in a cell, fig.3d), which is filled by a fluid. We used 
either air or a water-glycerol mixture at 60% concentration. The system is a harmonic oscillator with resonant 
frequency f a = \JC j I c g /(2tt) = oJq/ (27r) and a relaxation time r a = 2J e g jv = I /a. I c g is the total moment 
of inertia of the displaced masses (i.e. the mirror and the mass of displaced fluid) [11]. The damping has two 
contributions : the viscous damping v of the surrounding fluid and the viscoelasticity of the brass wire. 

The angular displacement of the pendulum 9 is measured by a differential interferometer [12, 13, 14, 15] 
which uses the two laser beams reflected by the mirror fig. 3a). The measurement noise is two orders of mag- 
nitude smaller than thermal fluctuations of the pendulum. 9{t) is acquired with a resolution of 24 bits at a 
sampling rate of 8192 Hz, which is about 40 times f a . We drive the system out-of-equilibrium by forcing it 
with an external torque M by means of a small electric current J flowing in a coil glued behind the mirror 
(Fig. 3b). The coil is inside a static magnetic field. The displacements of the coil and therefore the angular 
displacements of the mirror are much smaller than the spatial scale of inhomogeneity of the magnetic field. So 
the torque is proportional to the injected current : M = A.J ; the slope A depends on the geometry of the 
system. The practical realization of the montage is shown in figs. 3c), 3d). In equilibrium the variance 59 2 
of the thermal fluctuations of can be obtained from equipartition, i.e. 59 = Whs T/C ~ 2 nrad for our 
pendulum, where T is the temperature of the surrounding fluid. 

2.2 The equation of motion 

The dynamics of the torsion pendulum can be assimilated to that of a harmonic oscillator damped by the 
viscoelasticity of the torsion wire and the viscosity of the surrounding fluid, whose equation of motion reads in 
the temporal domain 

I cS 9+ f G(t-t')9(t')dt' + C9 = M + i 1 , (1) 

J — oo 

where G is the memory kernel and 7] the thermal noise. In Fourier space (in the frequency range of our interest) 
this equation takes the simple form 

[-I cS uj 2 + C]9 = M+ ri , (2) 

where : denotes the Fourier transform and C = C + i[C" + ojC^] is the complex frequency-dependent elastic 
stiffness of the system. C" and Cjf are the viscoelastic and viscous components of the damping term. 
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Figure 4: Equilibrium: The pendulum inside a glycerol-water mixture with M = 0. a) Square root of the 
power spectral density of 9. o directly measured spectrum, black dotted line is the spectrum estimated from the 
measure of x an d using eq.3 The red dashed and dotted lines show the viscous and viscoelastic component of 
the damping respectively, b) Probability density function of 9. The continuous line is a Gaussian fit 



2.2.1 Equilibrium 

At equilibrium,/. e. M = 0, the Fluctuation Dissipation Theorem (FDT) gives a relation between the amplitude 
of the thermal angular fluctuations of the oscillator and its response function. The response function of the 
system x = 0/M = -p* can be measured by applying a torque with a white spectrum. When M = 0, the am- 
plitude of the thermal vibrations of the oscillator is related to its response function via the fluctuation-dissipation 
theorem (FDT). Therefore, the thermal fluctuation power spectral density (psd) of the torsion pendulum reads 
for positive frequencies 



<|«|')_^lmS-^ . (3) 

The brackets are ensemble averages. As an example, the spectrum of 9 measured in the glycerol-water mixture 
is shown in fig.4a). In this case the resonance frequency is f a = y/ C / I c s / (2tt) = ojq/{2tt) = 217 Hz and 
the relaxation time r a = 2I c r/v = 1/a = 9.5 ms The measured spectrum is compared with that obtained 
from eq.3 using the measured x- The viscoelastic component at low frequencies correspond to a constant 
C'{ ^ 0. Indeed if u> — > then from eq.3 (\9\ ) oc 1/w as seen in fig.4a). Instead if C± = then for 
uj — > from eq.3 the spectrum is constant as a function of w. It is important to stress that in the viscoelastic 
case the noise 7] is correlated and the process is not Markovian, whereas in the viscous case the process is 
Markovian. Thus by changing the quality of the fluid surrounding the pendulum one can tune the Markovian 
nature of the process. In the following we will consider only the experiment in the glycerol-water mixture 
where the viscoelastic contribution is visible only at very low frequencies and is therefore negligible. This 
allows a more precise comparison with theoretical predictions often obtained for Markovian processes. The 
probability density function (pdf) of 6, plotted in fig.4b), is a Gaussian. 



2.3 Non-equilibrium Steady State (NESS): Sinusoidal forcing 

We now consider a periodic forcing of amplitude M Q and frequency ojd, i-e. M(t) = M o s'm(i0dt) [14]-[17]. 
This is a very common kind of forcing which has been already studied in the case of the first order Langevin 
equation [18] and of the two level system [19] and in a different context for the second order Langevin equation 
[20]. Furthermore this is a very general case because using Fourier transform, any periodical forcing can be 
decomposed in a sum of sinusoidal forcing. We explain here the behavior of a single mode. Experiments have 
been performed at various M a and w^. We present here the results for a particular amplitude and frequency: 
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M a = 0.78 pN.m and iOd/(^) = 64 Hz. This torque is plotted in Fig. 5a. The mean of the response to 
this torque is sinusoidal, with the same frequency, as can be seen in Fig. 5b. The system is clearly in a non- 
equilibrium steady state (NESS). 




Figure 5: a) Sinusoidal driving torque applied to the oscillator, b) Response of the oscillator to this periodic 
forcing (gray line) ; the dark line represents the mean response (9(t)). 

The work done by the torque M(t) on a time r n = 2ir n/uoa is 

W n = W T=Tn = / M(t)—dt (4) 

As 6 fluctuates also W n is a fluctuating quantity whose probability density function (pdf) is plotted in fig. 6a) 
for various n. This plot has interesting features. Specifically, work fluctuations are Gaussian for all values of 
n and W T takes negative values as long as r n is not too large. The probability of having negative values of 
W T decreases when r n is increased. There is a finite probability of having negative values of the work, in 
other words the system may have an instantaneous negative entropy production rate although the average of 
the work < W n > is of course positive (< . > stands for ensemble average). In this specific example is 
< W n >= 0.04 n(kB T). We now consider the energy balance for the system. 




|ff T (k„T/s) AL V ( k i3 T ) 

Figure 6: Sinusoidal forcing, a) Pdf of W T ; n = 7 (o), n = 15 (□), n = 25 (o) and n = 50 (x). b) Pdf of AC/ T 



2.4 Energy balance 

As the fluid is rather viscous we will take into account only the standard viscosity that is C'{ = and C'^ = v- 
In such a case eq. 1 simplifies 

I cS ^ + u— + C9 = M+ v , (5) 
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where i] is the thermal noise, which in this case is delta-correlated in time: < r/(t) rj(t') >= 2 T v5{t — t'). 

When the system is driven out of equilibrium using a deterministic torque, it receives some work and a 
fraction of this energy is dissipated into the heat bath. Multiplying Eq. (5) by 6 and integrating between t{ and 
ti + r, one obtains a formulation of the first law of thermodynamics between the two states at time ti and ti + r 
(Eq. (6)). This formulation has been first proposed in ref.[21] and used in other theoretical and experimental 
works [22, 18]. The change in internal energy AU T of the oscillator over a time r, starting at a time ti, is written 
as: 

AU T = U(U + r) - U(U) = W T -Q T (6) 



where W T is the work done on the system over a time r : 



W T = jf M(t') — (t')dt' 



(7) 



and Q T is the heat dissipated by the system. The internal energy is the sum of the potential energy and the 
kinetic energy : 



U{t) 



1 



r^cff 



d9 
dt 



(t) 



+ \ce{tf 



(8) 



The heat transfer Q T is deduced from equation (6) ; it has two contributions 



Q 7 



U+T 



I 2 



dt' 



tl+T dO 



(9) 



The first term corresponds to the viscous dissipation and is always positive, whereas the second term can be 
interpreted as the work of the thermal noise which has a fluctuating sign. The second law of thermodynamics 
imposes (Q T ) to be positive. 




Figure 7: Sinusoidal forcing, a) Pdf of W T b) Pdf of Q T for various n: n = 7 (o), n = 15 (□), n = 25 (o) and 
n = 50 (x). The continuous lines in this figures are not fits but are analytical predictions obtained from the 
Lnagevin dynamics as discussed in sect.3.4 



2.5 Heat fluctuations 

The dissipated heat Q T can not be directly measured because we have seen that eq.9 contains the work of the 
noise (the heat bath) that experimentally is impossible to measure, because rj is unknown. However Q T can 
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be obtained indirectly from the measure of W T and AU T , whose pdf measured during the periodic forcing are 
exponential for any r, as shown in fig.6b. We first do some comments on the average values. The average of 
AU T is obviously vanishing because the time r is a multiple of the period of the forcing. Therefore (W n ) and 
(Q n ) are equal. 

We rescale the work W T (the heat Q T ) by the average work (W T ) (the average heat (Q T )) and define: 
w T = -j^Tj (q T = 7qt)- m the present article, x T , respectively X T , stands for either w T or q T , respectively W T 
or Q T . 

We compare now the pdf of w T and q T in Fig. 7. The pdfs of heat fluctuations q n have exponential tails 
(Fig. 7b). It is interesting to stress that although the two variables W T and Q T have the same mean values 
they have a very different pdf. The pdf of w T are gaussian whereas those of q T are exponential. On a first 
approximation the pdf of q T are the convolution of a Gaussian (the pdf of W T ) and exponential (the pdf of 
AU T ). In Figs. 7 the continuous lines are analytical predictions obtained from the Langevin dynamics with no 
adjustable parameter (see sect.3.4). 



3 Fluctuation theorem 



In the previous section we have seen that both W T and Q T present negative values, i.e. the second law is 
verified only on average but the entropy production can have instantaneously negative values. The probabilities 
of getting positive and negative entropy production are quantitatively related in non-equilibrium systems by the 
Fluctuation Theorem (FTs). 

There are two classes of FTs. The Stationary State Fluctuation Theorem (SSFT) considers a non-equilibrium 
steady state. The Transient Fluctuation Theorem (TFT) describes transient non-equilibrium states where r mea- 
sures the time since the system left the equilibrium state. A Fluctuation Relation (FR) examines the symmetry 
around of the probability density function (pdf) p(x T ) of a quantity x T , as defined in the previous section. It 
compares the probability to have a positive event (x T = +x) versus the probability to have a negative event 
(x T = —x). We quantify the FT using a function S (symmetry function) : 

at x k B T fp(x T = +x)\ 

The Transient Fluctuation Theorem (TFT) states that the symmetry function is linear with x T for any values 
of the time integration r and the proportionality coefficient is equal to 1 for any value of r. 

S(x T ) = x T , \/x T , Vr. (11) 

Contrary to TFT, the Stationary State Fluctuation Theorem (SSFT) holds only in the limit of infinite time (r). 

lim S(x T ) = x T . (12) 

In the following we will assume linearity at finite time r [2, 29] and use the following general expression : 

S(x T ) = E x (r) x T (13) 

where for SSFT T, x (t) takes into account the finite time corrections and lim T _ i , 00 S x (r) = 1 whereas T, x (r) = 
1, Vr for TFT. 

However these claims are not universal because they depend on the kind of x T which is used. Specifically we 
will see in the next sections that the results are not exactly the same if X T is replaced by any one of W T , Q T and 
(T s to t :T ), defined in sect.3.3. Furthermore the definitions given in this section are appropriate for stochastic 
systems and in sect.6 we will discuss the differences between stochastic and chaotic systems. 
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3.1 Short history of FTs 



The first numerical evidence of relations of this kind has been given by Evans et al. in ref.[23] whereas the TFT 
was proved in ref.[24]. In 1995 Gallavotti and Cohen [25] proved SSFT for dynamical systems although in 
such a case x T takes a different meaning that we will discuss in sect.6. The proof of SSFT has been extended to 
stochastic dynamics in ref. [2, 26, 27, 28, 29]. Furthermore van Zon and Cohen proved that there is an important 
difference between the FTs for the injected power and those for the dissipated power [2, 29]. The SSFT has 
been proved also for other quantities such as the dissipation function [30] and the total entropy [31, 32]. Other 
theoretical papers studied FT and the reader may find a review in ref. [33, 34]. Experiments searching for FTs 
have been performed in dynamical systems [4, 5, 36], but interpretations are very difficult because a quantitative 
comparison with theoretical prediction can be doubtful. Other experiments have been performed in stochastic 
systems described by a first order Langevin equation: a Brownian particle in a moving optical trap [37] and 
an out-of-equilibrium electrical circuit [3] in which existing theoretical predictions [2, 29] were verified. Other 
experimental tests for FTs have been performed on driven two level systems [19] and on colloids [18]. 



3.2 FTs for Gaussian variables 

Let us suppose the the variable X T has a Gaussian distribution of mean < X T > and variance a\ r . It is easy 
to show that in order to satisfy FTs, the variable X T must have the following statistical property: 

a\ T =2k B T <X T > (14) 

This is an interesting relation because it imposes that the relative fluctuations of X T are 



a Xr 2k B T 



(15) 



x T V < x T > 

This means that the probability of having negative events reduces by increasing X T , specifically from eq. 15 
it follows that P(X T < 0) = erfc( y < X r > /(2 keT)) where erfc is the complementary error function. It 
is now possible to estimate the length of the time interval t i> s needed to observe at least one negative event, 
which is: 

r 

tobs = 7 — , s (16) 

where we used the fact that all the values W T computed on different intervals of length r are independent, 
which is certainly true if t is larger than the correlation time. 

Let us consider the specific example of section 2.3, i.e. < W T >= 0.04n(/csT) at u!d/{2 vr) = 67 Hz, 
M a = 0.78 pN.m and r = 2irn /oj^. The pdf of W T are Gaussian in this case (Fig. 6) and, as we will see in the 
next section, they satisfy SSFT for large r. Therefore eq.15 holds for X T = W T and we may estimate t b s in 
the asymptotic limit r Q < < r. For example at n = 200, one obtains from the above mentioned experimental 
values r ~ 3s >> r a and < W T >= 8ksT. Inserting these experimental values in eq.15 one gets roughly 
a negative event over an observational time t b s ~ 641s, which is already a rather long time for the distance 
between two events. For larger n and larger Mo this time becomes exponentially large. This justifies the fact 
that millions of data are necessary in order to have a reliable measure of SSFT. 



3.3 FTs for W T and Q T measured in the harmonic oscillator 

The questions we ask are whether for finite time FTs are satisfied for either x T = w T or x T = q T and what are 
the finite time corrections. In a first time, we test the correction to the proportionality between the symmetry 
function S(x T ) and x T . In the region where the symmetry function is linear with x T , we define the slope S x (r), 
i.e. S(x T ) = S x (r)x r . In a second time we measure finite time corrections to the value S x (r) = 1 which is 
the asymptotic value expected from FTs. 
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W T /<W T > q T Q V /<Q T > 

Figure 8: Sinusoidal forcing. Symmetry functions for SSFT. a) Symmetry functions S(w T ) plotted as a function 
of w T for various n: n = 7 (o), n = 15 (□), n = 25 (o) and n = 50 (x). For all n the dependence of S(w T ) 
on uv is lineal - , with slope X™(t). b) Symmetry functions S(q T ) plotted as a function of q T for various n. The 
dependence of S(q T ) on gr T is linear only for q T < 1. Continuous lines are is theoretical predictions. 
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Figure 9: Finite time corrections for SSFT. a) Sinusoidal forcing. T, w (t) as a function of n obtained from 
the slopes of the straight lines of Fig. 8a) (□). The circles correspond to another measurement performed at a 
different frequency. The finite time corrections depend on the driving frequency. The slope S g (r) measured for 
q T < 1 (Fig. 8b) have exactly the same values of S w (r) as a function of n. b) Linear forcing. E m (r) measured 
as a function of r with the driving torque M has a linear dependence on time. The finite time corrections 
depend on form of the driving. 

In this review article we will focus on the SSFT applied to the experimental results of sect. 2.3 and to other 
examples. The TFT will be not discussed here and the interested readers may look at ref.[17]. 

From the pdfs of w T and q T plotted in Figs.7, we compute the symmetry functions defined in eq. 10. The 
symmetry function S{w n ) are plotted in Fig. 8a) as a function of w n . They are linear in w n . The slope T, w (n) 
is not equal to 1 for all n but there is a correction at finite time (Fig. 9a). Nevertheless, T, w (n) tends to 1 for 
large n. Thus SSFT is satisfied for W T and for a sinusoidal forcing. The convergence is very slow and we have 
to wait a large number of periods of forcing for the slope to be 1 (after 30 periods, the slope is still 0.9). This 
behavior is independent of the amplitude of the forcing M Q and consequently of the mean value of the work 
(W n ), which, as explained in sec.3.2, changes only the time needed to observe a negative event. The system 
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satisfies the SSFT for all forcing frequencies uid but finite time corrections depend on ui^, as can be seen in 
Fig. 9a). 

We now analyze the pdf of q T (Fig.7b)) and we compute the symmetry functions S(q n ) of q n plotted in 
Fig. 8b) for different values of n. They are clearly very different from those of w n plotted in Fig. 8a). For 
S(q n ) three different regions appear: 

(I) For large fluctuations q n , S(q n ) equals 2. When r tends to infinity, this region spans from q n = 3 to 
infinity. 

(II) For small fluctuations q n , S(q n ) is a linear function of q n . We then define S 9 (n) as the slope of the 
function S(q n ), i.e. S(q n ) = S 9 (n) q n . We have measured [17] that S 9 (n) = E TO (n) for all the values of n, i.e. 
finite time corrections are the same for heat and work. Thus S g (n) tends to 1 when r is increased and SSFT 
holds in this region II which spans from q n = up to q n = 1 for large r. This effect has been discussed for the 
first time in refs.[29, 2]. 

(III) A smooth connection between the two behaviors. 

These regions define the Fluctuation Relation from the heat dissipated by the oscillator. The limit for large 
r of the symmetry function S(q T ) is rather delicate and it has been discussed in ref.[17]. 

The conclusions of this experimental analysis is that SSFT holds for work for any value of w T whereas for 
heat it holds only for q T < 1. The finite time correction to FTs, described by 1 — £ are not universal. They 
are the same both for w T and q T but they depend on the driving frequency as shown in Fig.9a). Furthermore 
they depend on the kind of driving force. In Fig. 9b) we plot T, w (r) measured when the harmonic oscillator is 
driven out of equilibrium by a linear ramp 1 . The difference with respect Fig. 9a) is quite evident. 

3.4 Comparison with theory 

This experimental analysis allows a very precise comparison with theoretical predictions using the Langevin 
equation (eq.5) and using two experimental observations: a) the properties of heat bath are not modified by 
the driving and b) the fluctuations of the W T are Gaussian (see also [38], where it is shown that in Langevin 
dynamics W T has a Gaussian distribution for any kind of deterministic driving force if the properties of the 
bath are not modified by the driving and the potential is harmonic). The observation in point a) is extremely 
important because it is always assumed to be true in all the theoretical analysis. In ref.[17] this point has been 
precisely checked. Using these experimental observations one can compute the pdf of q T and the finite time 
corrections S(r) to SSFT (see ref.[17]). The continuous lines in Fig. 9, Fig. 8 and Fig. 6 are not fit but analytical 
predictions, with no adjustable parameters, derived from the Langevin dynamics of eq.5 (see ref.[17] for more 
details). 

3.5 The trajectory dependent entropy 

In previous sections we have studied the energy W T injected into the system in the time r and the energy 
dissipated towards the heat bath Q T . These two quantities and the internal energy are related by the first law of 
thermodynamics (eq.9). Following notations of ref [31], we define the entropy variation in the system during a 
time r as : 



For thermostated systems, entropy change in medium behaves like the dissipated heat. The non-equilibrium 
Gibbs entropy is : 



where X t denotes the set of control parameters at time t and p(x(t),t, At) is the probability density function 
to find the particle at a position x(t) at time t, for the state corresponding to \ t . This expression allows the 

'The stationarity in the case of a ramp is discussed in ref.[29, 17] 
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Figure 10: a) Schematic diagram illustrating the method to compute the trajectory dependent entropy, b) Pdf of 
9{i) around the mean trajectory < 9(t) >. The continuous line is the equilibrium distribution 



definition of a "trajectory-dependent" entropy : 

s(t) = -k B \np(x(t),t,X t ) (19) 

The variation As to t,T of the total entropy s to t during a time r is the sum of the entropy change in the system 
during r and the variation of the "trajectory-dependent" entropy in a time r, As T = s(t + r) — s(t) : 

As to t,r = s to t(t + r) - Stot(*) = A%,r + As r (20) 

In this section, we study fluctuations of As totiT - computed using (17) and (19). We will show that As to t,r 
satisfies a SSFT for all r. In ref. [32], the relevance of boundary terms like As T has been pointed out for 
Markovian processes. 

We investigate the data of the harmonic oscillator described in sect.2.3.The probability to compute is the 
joint probability p(9(ti + T n ),9(ti + r n ), (p), where ip is the starting phase ip = tiUJd- The system is linear, so 
9(ti + T n ), B(ti + T n ) are independent; thus the joint probability can be factorized into a product. The expression 
of the trajectory dependent entropy is : 

As rn = - kB In ( Pm + r n )^ P (0(t l + r n ,p))\ ^ 

\p(9(ti + T n ),<p) p(0(ti + T n ,Lp) J 

For computing correctly the trajectory dependent entropy, we have to calculate the p(9(ti),ip) and p(9(ti),ip) 
for each initial phase <p (see Fig. 10a). These distributions turn out to be independent of ip and they correspond to 
the equilibrium fluctuations of 9 and 9 around the mean trajectory defined by (9(t)} and (9(t)). The distribution 
of 9{ti) is plotted in Fig. 10b), where the continuous line corresponds to the equilibrium distribution. Once the 
p(9(ti, ip) and p(9(ti), ip) are determined we compute the "trajectory-dependent" entropy. As fluctuations of 9 
and 9 are independent of tp we can average As Tn over <p which improves a lot the statistical accuracy. We stress 
that it is not equivalent to calculate first the pdfs over all values of ip — which would correspond here to the 
convolution of the pdf of the fluctuations with the pdf of a periodic signal — and then compute the trajectory 
dependent entropy. The results are shown in Fig. 11. 

In Fig. 11a), we recall the main results for the dissipated heat Q T = TAs m Tn . Its average value (T. As mjTn ) 
is linear in r n and equal to the injected work. The pdfs of T.As m Tn are not Gaussian and extreme events 
have an exponential distribution. The pdf of the "trajectory-dependent" entropy is plotted in fig. lib); it is 
exponential and independent of n. We supeipose to it the pdf of the variation of internal energy divided by T at 
equilibrium: the two curves match perfectly within experimental errors, so the "trajectory-dependent" entropy 
can be considered as the entropy exchanged with the thermostat if the system is at equilibrium. The average 
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Figure 11: Torsion pendulum, a) pdfs of the normalized entropy variation As m]Tn /(As mjTn ) integrated over 
n periods of forcing, with n = 7 (o), n = 15 (□), n = 25 (o) and n = 50 (x). b) pdfs of As Tn , the 
distribution is independent of n and here n = 7. Continuous line is the theoretical prediction for equilibrium 
entropy exchanged with thermal bath As m TniCq . c) pdfs of the normalized total entropy As t ot,r n /(Astot,r n }> 
with n = 7 (o), 7?, = 15 (□), n = 25 (o) and n = 50 (x). d) Symmetry functions for the normalized entropy 
variation in the system (small symbols in light colors and X T stands for T As m>Tn = Q T ) and for the normalized 
total entropy (large symbols in dark colors and X T stands for T As to t,r n ) for the same values of n. 



value of As Tn is zero, so the average value of the total entropy is equal to the average of injected power divided 
by T. In Fig. 11c), we plot the pdfs of the normalized total entropy for four typical values of integration time. 
We find that the pdfs are Gaussian for any time. 

The symmetry functions (eq. 10) of the dissipated heat S(TAs m )Tn = Q T ) and the total entropy 5(TAs totiTn ) 
are plotted in Fig. 1 Id). As we have already seen in Fig.8, S(Q T ) is a non linear function of Q T = T As m)T . 
The linear behavior, with a slope that tends to 1 for large time, is observed only for for As m Tn < (As m Tn ) < 1. 
For the normalized total entropy, the symmetry functions are linear with As t ot,r n f° r ai l values of As t ot,T n an d 
the slope is equal to 1 for all values of r n . Note that it is not exactly the case for the first values of r n because 
these are the times over which the statistical errors are the largest and the error in the slope is large. 

For the harmonic oscillator we have obtained that the "trajectory-dependent" entropy can be considered as 
the entropy variation in the system in a time r that one would have if the system was at equilibrium. Therefore 
the total entropy is the additional entropy due to the presence of the external forcing : this is the part of entropy 
which is created by the non- equilibrium stationary process. The total entropy (or excess entropy) satisfies the 
Fluctuation Theorem for all times and for all kinds of stationary external torque[31, 32]. More details on this 
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problem can be found in ref.[39]. 



4 The non-linear case: stochastic resonance 




Figure 12: a) Drawing of the polystyrene particle trapped by two laser beams whose axis distance is about the 
radius of the bead, b) Potential felt by the bead trapped by the two laser beams. The barrier height between the 
two wells is about 2k B T. 

The harmonic oscillator cannot be driven to a non linear regime without forcing it to such a high level where 
thermal fluctuations become negligible. Thus in order to study the non linear effects we change experiment 
and we measure the fluctuations of a Brownian particle trapped in a non-linear potential produced by two 
laser beams, as shown in Fig. 12. It is very well known that a particle of small radius R ~ 2 fim is trapped 
in the focus of a strongly focused laser beam, which produces a harmonic potential for the particle, whose 
Brownian motion is confined inside this potential well. When two laser beams are focused at a distance D ~ R, 
as shown in Fig. 12a) the particle has two equilibrium positions, i.e. the foci of the two beams. Thermal 
fluctuations may force the particle to move from one to the other. The particle feels an equilibrium potential 
Uq{x) = ax 4 — bx 2 — dx, shown in Fig. 12b), where a, b and d are determined by the laser intensity and 
by the distance of the two focal points. This potential has been computed from the measured equilibrium 
distribution of the particle P(x) oc exp(Uo(x)). The right left asymmetry of the potential (Fig. 12b) is induced 
by small unavoidable asymmetries, induced by the optics focusing the two laser beams. In our experiment the 
distance between the two spots is 1.45 /mi, which produces a trap whose minima are at x m i n = ±0.45 /im. 
The total intensity of the laser is 29 mW on the focal plane which corresponds to an inter-well barrier energy 
5U a = 1.8 k B T, ax 4 mm = 1.8 k B T, bx 2 mm = 3.6 k B T and d\x min \ = 0.44 k B T (see ref.[40] for more 
experimental details). The rate at which the particle jumps from one potentials well to the other is determined 
by the Kramer's rate = ^- exp( ) where r Q is a characteristic time. In our experiment ~ 0.3 Hz at 
300 K. 

To drive the system out of equilibrium we periodically modulate the intensity of the two beams at low 
frequency. Thus the potential felt by the bead is the following profile: 

U(x,t) = U Q {x) + U p (x,t) = U + cx sin(27r/i), (22) 

with c\x m i n \ = 0.81 k B T. The amplitude of the time dependent perturbation is synchronously acquired with 
the bead trajectory. 2 

2 The parameters given here are average parameters since the coefficients a, b and c, obtained from fitted steady distributions at 
given phases, vary with the phase (8a/ a ~ f 0%, Sb/b w 8c/ c ~ 5%). 
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Figure 13: a) The perturbed potential att = jj and t = jj. b) Example of trajectory of the glass bead and 
the corresponding perturbation at / = 0.1 Hz. c) Mean injected energy in the system over a single period 
as a function of the driving frequency. < W s > □ and < W c i > o coincide as their mean values are equal 
within experimental errors. The error bars are computed from the standard deviation of the mean over different 
runs. Inset: Standard deviations of work distributions over a single period normalized by the average work as a 
function of the frequency (same symbols). 



An example of the measured potential for t = jj and jj is shown on the Fig. 13a). This figure is obtained 
by measuring the probability distribution function P(x,t) of x for fixed values of csm(2irft), it follows that 
U(x,t) = -]n(P(x,t)). 

The x position of the particle can be described by a Langevin equation: 

dU(x,t) 

ix = — h tj, (23) 

with 7 = 1.61 10~ 8 N s m _1 the friction coefficient and r\ the thermal noise delta correlated in in time. When 
c ^ the particle can experience a stochastic resonance [41], when the forcing frequency is close to the 
Kramer's rate. An example of the sinusoidal force with the corresponding position are shown on the figure 
13b). Since the synchronization is not perfect, sometimes the particle receives energy from the perturbation, 
sometimes the bead moves against the perturbation leading to a negative work on the system. Two kinds of 
work can be defined in this experiments [40] 



^tHjT'*^ ( 24 ) 

W d , n (t) = - J™' dtx^^ (25) 

(26) 

where in this case tf = j is a multiple of the forcing period. The work W S)n is the stochastic work (used in 
Jarzynsky and Crooks relations [47, 49, 13]) and Wd, n is the classical work that will be discussed in this article. 
The results on W Sin are quite similar but there are subtle differences discussed in ref.[40]. 

We first measure the average work received over one period for different frequencies (tf = 4 in eq. 25). 
Each trajectory is here recorded during 3200 s in different consecutive runs, which corresponds to 160 up to 
6400 forcing periods, for the range of frequencies explored. In order to increase the statistics we consider 
10 5 different t D . The figure 13c) shows the evolution of the mean work per period for both definitions of the 
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Figure 14: a) Distribution of classical work W c i for different numbers of period n = 1, 2, 4, 8 and 12 (/ = 
0.25 Hz). Inset: Same data in lin-log. b) Normalized symmetry function as function of the normalized work 
forra = 1 (+), 2 (o), 4 (o), 8 (A), 12 (□). 



work. First, the input average work decreases to zero when the frequency tends to zero. Indeed, the bead hops 
randomly several times between the two wells during the period. Second, in the limit of high frequencies, the 
particle has not the time to jump on the other side of the trap but rather stays in the same well during the period, 
thus the input energy is again decreasing when increasing frequency. In the intermediate regime, the particle 
can almost synchronize with the periodical force and follows the evolution of the potential. The maximum of 
injected work is found around the frequency / 0.1 Hz, which is comparable with half of the Kramers' rate 
of the fixed potential tk = 0.3 Hz. This maximum of transferred energy shows that the stochastic resonance 
for a Brownian particle is a bona fide resonance, as it was previously shown in experiments using resident time 
distributions [42, 43] or directly in simulations [44, 45]. In the inset of Fig. 13, we plot the normalized standard 
deviation of work distributions (a/ (W)) as a function of the forcing frequency. The curves present a minimum 
at the same frequency of 0. 1 Hz, in agreement again with the resonance phenomena. 

In order to study FT for stochastic resonance we choose for the external driving a frequency / = 0.25 Hz, 
which ensures a good statistic, by allowing the observation of the system over a sufficient number of periods. 
We compute the works and the dissipation using 1.5 10 6 different t on time series which spans about 7500 
period of the driving. 

We consider the pdf P(W c i) which is plotted in (Fig. 14a). Notice that for small n the distributions are 
double peaked and very complex. They tend to a gaussian for large n (inset of Fig. 14a). On Fig. 14b), we 
plot the normalized symmetry function of W c i, n . We can see that the curves are close to the line of slope one. 
For high values of work, the dispersion of the data increases due to the lack of events. The slope tends toward 
1 as expected by the SSFT It is remarkable that straight lines are obtained even for n close to 1, where the 
distribution presents a very complex and unusual shape (Fig. 14a). We do not discuss here the case of W SjTl as 
the behavior is quite similar to that of Wd,n [40]. The very fast convergency to the asymptotic value of the the 
SSFT is quite striking in this example. The measurement are in fully agreement with a realistic model based on 
the Fokker Planck equations where the measured values of U(x, t) has been inserted [46]. This example shows 
the application of FT in a non-linear case where the distributions are strongly non-Gaussian. 

5 Applications of Fluctuation Theorems 

The Fluctuations Theorems have several important consequences such as the Jarzinsky and Crooks equalities[47, 
48, 49], which are useful to compute the free energy difference between two equilibrium states using any kind 
of transformation[13, 15, 51, 50]. The Hatano-Sasa[22] relations and the recently derived Fluctaution Dissipa- 
tion Theorems[52] are related to FTs and are useful to compute the response of a NESS using the steady state 
fluctuations of the NESS. As we have seen the FT allows the calculation of tiny amount heat, which can be 
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useful in many applications in aging systems [53, 54] and biological systems. 

The FTs for Langevin systems can be used to measure an unknown averaged power. This idea has been 
discussed first in the context of electrical circuits [3] and in ref.[57] it has been applied for the first time to the 
measure of the the torque of a molecular motor. We discuss the method in some details in the next subsection. 
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Figure 15: Molecular motor, a) Schematic diagram of a Fx-ATPase molecular motor composed by a rotor 7 
(radius ~ 1 nm) which rotates inside stator of radius ~ 5 nm formed by three a /3 subcomplexes. Sequential 
chemical reactions between the stator and the rotor produce the motion. The a/3 subcomplexes are attached to 
a suitably coated glass plate. Streptavidin is used to attach to 7 either actin filaments [58] or streptavidin-coated 
beads [57]. b) In order to follow the rotation of the rotor with a standard microscope a streptavidin-coated bead 
of radius ~ 0.5 /mi is glued to 7 subunit (drawing not to scale). The figures c),d) and e) (taken from ref.[57])) 
illustrate the results of a measure, c) Pdfs for several r of A9 T . d) Symmetry function extracted from the Pdfs 
of c). e) Slopes NTj(t) of the symmetry function as a function of r. The different colors pertain to different 
experimental conditions. Notice the convergency to an unique value of N for large r. Strictly speaking in this 
figure the function S(r) keeps into account also the fact that eq.27 is not necessarily valid for short times (see 
text) 



5.1 Measuring the power of a molecular motor 

A molecular rotary motor, as any kind of motor, is constituted by a stator and a rotor. The movement of the 
rotor is provoked by chemical reactions occurring sequentially between the rotor and stator. A typical example 
of bio-motor is the bacterial flagellum. However it has been shown [58] that a single molecule of Fi-ATPase 
may act as a motor composed by a 7 subunit ( radius ~ 1 nm) which rotates inside a barrel of radius ~ 5 nm 
formed by three a /3 subcomplexes (see fig. 15a) for a schematic diagram and ref.[57, 58] for more details). In 
experiments, the a(3 subunit is stuck on a suitably activated glass plate as shown in fig. 15a). The measure of 
the torque of this motor is important in order to know its efficiency as a function of the concentration of the 
chemicals contained in the liquid surrounding it. The typical size of this molecular motor is several nanometers 
and the moving unit is too small to be observed with an optical microscope. Therefore to measure the torque of 
Fi-ATPase motor a streptavidin-coated bead of radius 0.5 fim is glued on the subunit 7, and the motion of this 
bead is followed by a standard microscope as sketched in Fig. 15b). The motion of the bead occurs on a torus 
and 9 is the coordinate of the motion along the torus. The time evolution of 9 can be described by a Langevin 
equation: 
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T6 = N + t 



(27) 



where N is the torque of the motor, F the viscous dissipation of the fluid surrounding the bead and £ the 
thermal noise. In order to estimate N, it is in principle enough to measure the mean angular velocity < 9 > of 
the bead and obviously from eq.27 one gets: 

N = T <9> (28) 

However the procedure is not so simple because in order to estimate N one has to know the exact value of T, 
which is a function of the viscosity of the fluid, the radius L of the bead path, the radius R of the bead, the 
shape of the bead and the distance Z Q of the bead from the surface of the glass plate where the experiment 
is performed (see Fig. 15b). The variable Z Q is certainly the most difficult to be determined experimentally. 
Therefore using eq.28 the error on N can be really very large. There is instead another method proposed in 
ref.[3] for electric circuits and first used for measurements in bio-motors in ref.[57]. This method is based 
on work fluctuations and is much more precise. To apply it, we suppose that N is constant. This is a quite 
reasonable hypothesis for certain regimes of Fi-ATPase motors. We compute the work W T performed by the 
motor in a time r: 

rt+r 

W T {t) = J N 0(t)dt = N A9 T (29) 

where A9 T = 9(t + r) — 9(t) and we have used the fact that N is constant. In previous section we have seen 
that Langevin systems satisfy the SSFT, which we now apply to W T . Using eq.29 and the constancy of N SSFT 
for the molecular motor reads: 

In ( P ,l\ ) = E(t) N — with S(r) -> 1 for r -> oo (30) 

\P(-A9 T ) J y ' k B T y 1 

This equation is quite interesting because the value of N can be determined only by the measure of the 
fluctuations of A9 T . Indeed plotting ln(P(A0 r )/P(— A0 T )) as a function of A9 T /(k B T) we notice that the 
slope of the straight lines is S(r) N. Therefore studying the asymptotic value of this slope for large r one can 
determine N. It is interesting to note that in this case the knowledge of T is not needed. This technique has 
been recently applied to molecular motor in ref.[57] and their main results are plotted in Figs. 15 c),d),e). The 
relevant parameter S(r) N, extracted from the pdf of A9 T (Fig. 15 c) and the symmetry function (Fig. 15 d), 
is plotted in Fig. 15e) 3 . We see a clear convergence to a unique value and one gets a very precise estimation 
of the torque of the molecular motor independently of the size and shape of the bead glued to 7 unit of the 
Fi-ATPase motor. This is a very specific and interesting example of the possible applications of FT. 

6 The chaotic systems 

In previous sections we have studied the probability of the instantaneous negative entropy production rates 
within the context of the FTs for stochastic systems, where the fluctuations are produced by the coupling with 
a thermal bath. In sec. 3.2 we have seen that when the energy injected into the system is larger than 100 k B T 
the probability of these negative events is very small and the time needed to observe them becomes extremely 
long. In other words the role of thermal fluctuations becomes negligible. 

However in the introduction we have shown that instantaneous negative entropy production rates can be 
observed in chaotic systems such as, for example, turbulence and granular media, where the fluctuations are 
produced by the non linear interactions of many degrees of freedom. We have also pointed out that for chaotic 
systems the amount of injected energy is order of magnitudes larger than k B T and of course thermal fluctuations 

3 It has to be stressed that in this specific case E(r) keeps into account the fact that for short time eq.27 is not a good approximation 
for the dynamics of the motor and eqs.29,30 apply only for long time 
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do not play any role in the fluctuating dynamics. The question that we want to analyze in this section is whether 
we can apply in these systems the FTs defined in sec. 3 for stochastic systems, eqs. 10-13. For a dissipative 
chaotic system one could imagine to replace &bT, in eqs. 10-13, with a characteristic energy E c which keeps 
into account the relevant energy scales of the system fluctuations. However the definition of this relevant energy 
scale can be in general difficult and even impossible, because it may depend on the observable and on the kind 
of forcing. Thus the approach of introducing an E c is not very useful to compare the experimental results with 
the proof given for dynamical systems [25]. Indeed in this case the theorem considers a quantity : 

Or l! +T o(t) dt 

y T = —2— = h ^ (31) 

where a(t) is the instantaneous phase space contraction rate, a T the integral of a on a time r and < a T > the 
mean of a T . Three hypothesis has been done on the dynamical system which must be : a) dissipative, b) time 
reversible c) Anosov 4 

1 P(v ) 

- In — — — - = < a > y T + Ofl/r) for r -> oo (32) 

where < a > is the mean phase space contraction rate, which has the dimension of 1/t. In this equation 
the relevant variable is the phase space contraction rate which has been identified as the entropy production 
rate[25]. The phase space contraction rate is a global variable of the system but an extension of the theorem for 
local variables has been done in refs.[55, 56]. One reason for developing local FT is that global fluctuations are 
usually not observable in macroscopic systems, as a consequence it is important to understand whether a local 
measurement is representative of the dynamics. Eq.32 has been tested in several numerical simulations (see for 
example ref.[33, 34] for a review), here we want to focus on experiments. 



6.1 Experimental test 

The test of eq.32 in experiments is extremely useful to analyze several important questions. The first one is 
whether eq.32 may have a more general validity independently of the restrictive hypothesis done to prove it. 
Indeed the hypothesis b) is never satisfied in real systems and the hypothesis c) does not necessarily apply to all 
of them. Thus in general we do not even know whether eq.32 can be applied in the experimental system under 
study. The second question concerns the choice of the observable. Indeed the direct measure of the phase space 
contraction rate is not possible and one has to rely upon the measure of another observable usually the energy 
W T injected into the system by the external forces in a time r. In other words one is making the important 
hypothesis that y T , defined in eq.31 in terms of a(t) is equivalent to x T = W T / < W T >. This hypothesis, that 
is not necessarily valid, is the second question that one would like to address in experiments. The third question 
is related to the estimation of the prefactor < a > in the right hand side of eq.32. This prefactor, which is a 
function of the Lyapunov exponents, is very difficult to estimate in an experimental system. Finally the last 
question concerns the relevance of a local observable to characterize the dynamics of the system. 

There are not many experiments where these questions have been analyzed in some details. In several 
experiments [4, 36] only the linearity in x T of the symmetry function p{x T ) = (1/t) In (P(x T ) / P(—x T )) has 
been checked, which, for the reasons discussed in the previous paragraph, is only a partial test. For example 
three experiments have tried to give an answer to the question of the prefactor. Two of these experiments are 
performed in granular media [5, 6] and the third on mechanical wave turbulence in a metallic plate[63]. We 
will not describe in details the experiments here but we will comment the main results. 



6.1.1 Granular media 

The two experiments of refs.[5, 6] consist of diluted granular media strongly shaken by a vibrator, but the 
measured quantities are not the same. In ref.[5] the authors measure the fluctuations of the energy flux in a 

4 For a precise definition see refs.[35, 34]. Roughly speaking this property ensures that the system is chaotic and that on the attractor 
there are no regions of finite volume that do not contain points. 



19 



subvolume of the system. Instead ref . [6] the work done by an external force on a ratchet inside the granular 
media is measured. They both find that although the system is not thermal the stochastic version of SSFT 
(eq.10) holds, provided that ksT is replaced by a characteristic energy E c : 



In 



PjXr) 
P(-Xr) 



<X T > 



x t + 0{1/t). 



(33) 



with x T = X T j < X T >, X T is the integral of energy flux in ref.[5] and the work W T performed by an 
external force on a ratchet in ref. [6]. It must be pointed out that in both experiments of refs.[5, 6] the energy 
E c has been measured independently. In ref. [5] is found that E c is about 5 times larger than the kinetic energy 
Ke of the shaken granular medium for all the values of the control parameters used in the experiment. The 
fact that Ke and E c have the same dependence on the control parameter have been interpreted considering 
that the vibrator injects into the system the amount of energy lost in the collisions but once excited in a NESS 
the granular medium behaves like a thermal bath for the measured observable. However this interpretation is 
not necessarily correct. Indeed for the experiment of ref. [5] it as been shown in a numerical simulation [59] 
that for the quantity measured in this experiment FT does not apply for the large deviations x T because the 
symmetry function S(x T ) becomes non-linear for large x T . This discrepancy between theory and experiment 
is obviously coming from the fact that experimentally the very large deviation are difficult to be measured, thus 
the non-linear part of S(x T ) cannot be observed. However the experiment of ref. [5] is certainly interesting 
because is the first where the question of the prefactor has been analyzed experimentally. In contrast for the 
experiment of ref. [6] is observed that E c = Ke{1 + o>)/2 where a is the restitution coefficient of the grains. 
It is interesting to notice that Ke(1 + a)/2 is the temperature of an intruder inside a diluted granular gas as 
it has been found in theoretical models [60, 61, 62]. These two examples of comparison between numerical 
and experimental results show the difficulty of interpreting the experimental results on FT and the importance 
of verifying them in a precise theoretical framework. Where does the difference between the experiments of 
ref. [5] and ref. [6] come from ? The answer can be found on the fact that in the two experiments two different 
quantities are measured. Indeed in sect.3 we have seen that, even for stochastic systems, the fluctuations of 
W T , Q T and As to t,T n behave differently within the context of FT. This can be more complex for granular media 
and it will be useful to give more insight on this point. Furthermore even in cases where a description in terms 
of Ec applies, comparing eq.33 with eq.32 a question that arises naturally is whether < W T > /E c is a good 
estimation of < a > r. This is an important question which will be interesting to analyze in the future. 

6.1.2 Mechanical waves 

In the experiment on mechanical waves [63], a metallic plate is set into a chaotic state of wave turbulence by a 
periodic local forcing at 75 Hz. (see Fig. 16a). The chaotic dynamics is produced by the non-linear interaction 
of the oscillatory modes of the plate. The authors measure the local force and displacement (see Fig. 16a) and 
compute the work W T done on a time r by the external force which excites the vibrations of the plate. They 
find that the pdfs of x T = W T / < W T > are strongly non-gaussian (see Fig. 16b). From these pdfs they 
compute the symmetry function p(x T ) which is plotted (see Fig. 16c) as a function of x T . We see that in spite 
of the fact that the pdf are not Gaussian the function p(x T ) (Fig. 16c) converges to a unique straight line for 
large r as predicted by FT. From eq.32 the slope of this straight line is < a >, which the authors can estimate 
independently by measuring the relaxation time of the vibrational modes. They find that the values estimated 
with the two methods (FT and the relaxation time) are very close and within experimental errors. This result is 
quite interesting and it is probably the only experiment where a direct test of eq.32 has been done. Certainly 
the errors of this comparison are very large but this kind of tests are useful to understand in some details the 
applications of FTs to chaotic systems. 
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Figure 16: Mechanical waves in a metallic plate (from ref.[63]). a) Schematic diagram of the experiment. 
A steel plate is suspended to the frame. Dimensions are in cm. The electromagnetic exciter produce et lo- 
cal forcing of the plate vibrations. A laser vibrometer measures the normal velocity at the excitation point. 
b)and c) Results on the periodic forcing, (b) Pdfs of the injected power on the time durations of r for 
r = 3.5ms, 6.5, 13.5, 20, 26.5, 33.5, 40, 47.5, 52ms ) (c) Functions p(x T ) = (1/r) In (P(x T )/P(-x T )) 
obtained from the Pdfs of (b). Inset: compensated value p(x T )/x T in a semi log plot. In (c) the dashed line 
corresponds to a linear law of slope < a >= 700Hz. 



7 Summary and concluding remarks 

In this paper we have reviewed several experimental results on the fluctuations of injected and dissipated power 
in out of equilibrium systems. We considered the two cases when the fluctuations are produced by the coupling 
with the heat bath (stochastic systems) and when they are produced by the non linear interactions of many 
degrees of freedom (chaotic systems). We have seen that in both cases we observe that the external forces may 
produce a negative work because of fluctuations. The probability of these negative events has been analyzed in 
the framework of fluctuation theorem. 

We have mainly discussed the stochastic systems described by Langevin equations, both with harmonic and 
unharmonic potential. We have seen that injected and dissipated power present different behaviors. FTs are 
valid for any value of W T whereas can be applied only for Q T < < Q T > in the case of the heat. We have 
also seen the the finite time corrections to SSFT depend on the driving and on the properties of the system. 
We have introduced the total entropy, which takes into account only the entropy produced by the external 
forces neglecting the the equilibrium fluctuations. For the total entropy FTs are valid for all the times. We 
discussed the applications of FTs to extract important physical properties of a stochastic system. Thus one may 
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conclude that for Markovian systems driven by a deterministic force the applications of FT does not present any 
major problems and can be safely applied. The case of random driving has been recently discussed and several 
problems may arise when the variance of the driving become larger than the fluctuations induced by the thermal 
bath. We have not discussed this problem but an analysis of this specific case can be found in refs.[8, 64, 65]. 

Finally we discussed the applications of FT to chaotic systems. The experimental test is in this case very im- 
portant and useful because many questions can be asked on the system under study which does not necessarily 
verify all the theoretical hypothesis. One has to say that in the case of non-Gaussian statistics even the linearity 
of the symmetry function can be an interesting result. However we pointed out that, for a real comparison with 
theory, the difficulty is to estimate of the prefactor of eq.32 by an independent measurement. Only a few exper- 
iments have addressed this point in some details, but many problems remain open and it seems to be difficult to 
find a general behavior for chaotic systems as for the case of stochastic ones. 
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